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1.  ABSTRACT 

Tools  for  electromagnetic  simulation,  such  as  SPECRAY  EM  /  FERMAT,  make  RCS  and  ISAR  image 
computations  down  to  millimetre  wavelength  possible  from  the  CAD  of  targets.  The  higher  the  frequency  is, 
the  greater  the  influence  of  the  precision  and  the  realism  of  the  CAD  models  on  results  are.  In  millimetre 
wavelength,  whatever  the  precision  of  the  CAD,  the  roughness  of  surface  and  the  flatness  of  sheets  of 
construction  of  terrestrial  military  vehicles  cannot  be  reproduced  by  the  geometry.  Nevertheless,  a  realistic 
prediction  of  RCS  can  be  carried  out  by  taking  into  account  the  statistical  surface  qualities,  which  are 
superposed,  on  the  deterministic  effects  of  the  geometry  of  the  CAD.  The  calculations  carried  out  on  the  CAD 
600  and  80  kfacets  of  the  ZSU  23-4  illustrate  these  effects,  without  degrading  the  efficiency  of  computations. 

2.  THE  ELECTROMAGNETIC  COMPUTER  CODE  FERMAT 

The  SPECRAY  EM  software  includes  several  tools  for  the  generation,  management  of  physically  textured 
databases,  the  edition  of  scenarios  relating  to  a  given  sensor,  the  post-processing  and  visualisation  of  the 
output  data  and  an  EM  core,  FERMAT.  FERMAT  deals  with  the  electromagnetic  calculations  of  the 
interactions  between  EM  wave  and  the  observed  scene. 

The  FERMAT  software  has  been  developed  by  the  Electromagnetism  and  Radar  Department  of  ONERA 
within  the  framework  of  a  project  dedicated  to  electromagnetic  (EM)  modelling.  This  computer  code 
capitalises  the  results  of  research  activities  in  high  frequency  scattering  methods.  A  unique  tool  is  currently 
available  for  studies  dealing  with  radar  analysis,  antennas  radiation,  EM  inter-system  compatibility  and 
propagation  applications.  This  numerical  simulation  tool  has  the  ambition  to  calculate  scattered  EM  fields  at 
high  frequencies  (i.e.  the  size  of  objects  is  supposed  large  compared  to  the  wavelength),  in  a  virtual  3  D, 
geometrical  and  physical  complex  environment  including  natural  and  man  made  objects.  EM  fields 
calculations  either  on  a  surface  or  in  a  volume  could  be  carried  out  according  to  specific  applications. 

FERMAT  associates  various  techniques  and  tools. 

Virtual  3D,  geometrical  data  bases,  composed  of  a  great  number  of  elements,  landscapes  and  objects  such 
as  vehicles,  buildings,  ...,  are  modelled  by  thousands  of  plane  geometrical  faces  and  associated  management 
tools.  Specific  features  and  textures  relative  to  electromagnetic  scattering  are  provided  for  each  face. 
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Figure  1:  Examples  of  natural  and  man  made  databases 


A  geometrical  Shooting  and  Bouncing  Rays  (SBR)  technique  has  been  optimised  in  order  to  calculate  the 
intersections  between  rays  from  the  transmitter  towards  the  database  and  back  to  a  receiving  point.  EM  models 
of  propagation,  reflection,  diffraction  and  an  operating  strategy  (thanks  to  SBR)  which  allows  unified 
calculation  for  the  near  or  far  EM  scattered  fields  from  the  scenes.  These  models  are  the  formulations  of 
Geometrical  Optics  (GO),  Physical  Optics  (PO),  Uniform  Theory  of  Diffraction  (UTD)  and  Equivalent 
Current  Method  (ECM).  Again,  these  models,  judiciously  coupled  with  the  SBR  technique  make  the 
computation  time  slightly  dependent  on  the  complexity  of  the  treated  virtual  scene. 


Figure  2:  Coupling  of  SBR  technique  and  EM  asymptotic  methods 


Generation  and  management  scenarios  tools  has  been  developed  to  simultaneously  simulate,  for  example, 
motions  of  a  radar  sensor  and  displacements  of  mobiles  in  the  landscape  according  to  the  type  of  field 
calculations.  Outputs  and  visualisation  tools  have  been  adapted  to  investigated  applications. 
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Figure  3:  Examples  of  ISAR  images  from  the  canonical  target  set  and  the  ZSU  23-4 


The  strength  of  FERMAT  is  that  it  is  developed,  under  an  agreement  of  partnership  between  ONERA  and 
OKTAL  SE.  This  company  brings  its  competencies  in  the  management  and  exploitation  of  databases,  and  the 
control  of  the  SBR  technique. 

The  validation  mainly  consisted  in  comparison  with  the  EM  models  calculated  either  with  other  reference 
methods  or  with  measurements  on  simple  object  known  as  canonical  and  allowing  evaluating  the  various  EM 
interactions.  These  comparisons  took  place  during  specialised  national  and  international  "workshops". 


3.  ROUGHNESS  CHARACTERISTICS 


3.1.  Introduction 

In  previous  investigation  related  to  millimetre  wavelength  modelling,  under-meshing  of  objects  had  led  to 
peaks  of  scattering  in  the  specular  direction  (very  high  values  RCS  that  do  not  correspond  to  predicted  or 
known  values  from  real  targets).  A  disturbance  of  height  on  the  object  facets  was  introduced  to  correct  this 
effect,  by  attenuating  it.  The  technique  consisted  in  plating  a  texture  (space  height  function)  of  random  height 
on  the  smooth  facets  of  the  under-meshed  target.  We  expected  that  the  technique  also  made  possible  to  take 
into  account  the  presence  of  real  roughness  and  non-flatness  of  real  surfaces  giving  results  superimposing 
random  effects  on  the  coherent  scattering. 

ONERA,  OKTAL  SE  and  CLEAR  completed  a  work  of  expertise,  in  millimetre  wavelength,  to  validate 
those  first  ideas  in  a  context  representative  of  the  requirements.  Measurements  on  a  VAB  vehicle  (made  on 
CLEAR  STRADI  set  up  in  Rennes,  France)  and  computation  (by  FERMAT)  on  a  CAD  file  of  this  same 
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vehicle  equipped  with  a  texture  of  roughness  have  been  compared.  The  principal  conclusions  were  a  better 
agreement  between  the  predictions  and  the  importance  of  the  presence  roughness  texture  on  the  surface  of  the 
object.  Through  this  expertise,  ONERA  proved  that  the  method  is  well  adapted  to  simulation  of  radar 
behaviour  radar  in  millimetre  wavelength. 


ISAR  images  at  millimetre  wavelength 

Angular  sector  -0.5/-0.5  deg 

Measurements  Calculations 


Figure  4:  Previous  results  in  using  roughness  texture  on  a  recovery  vehicle 


3.2.  Profile  parameters 

Generally,  a  surface  roughness  is  characterised  by  the  root  mean  square  a  of  its  height  distribution  relating 
to  wavelength  and  by  the  space  correlation  length  Lc  of  an  extracted  profile  from  this  surface.  The  interaction 
of  EM  waves  and  the  surface  presents  a  coherent  part  in  the  specular  direction  and  a  scattering  part  in  the 
other  directions.  According  to  appropriate  bibliography,  the  calculation  of  the  specular  component  does  not 
require  the  knowledge  of  correlation  length.  However,  this  one  is  necessary  an  evaluation  of  the  scattering 
component. 

A  surface  is  known  as  smooth  when  the  product  ka  <0.2  or  —<0.03. 

X 

A  surface  is  known  as  rough  when  the  product  k(j>l  or  —  >0.16  with  k  = - 

X  X 

According  to  Rayleigh  [Ulaby  and  A1  82],  the  coherent  component  is  present  for  —  <  0.25 . 

X 

Other  characterisation  criteria  using  the  correlation  length  will  be  specified  thereafter.  Kirchhoff 
approximation  [Ulaby  and  A1  82]  led  to  an  expression  of  the  attenuation  of  the  coherent  scattering  (specular 
reflection)  in  the  form  known  as  the  modified  Fresnel  reflection  coefficient: 

R  =  R  '  exp(-2k^  cos^  ( 6^ )) 

Note:  Thereafter,  we  will  use  the  term  of  "bump"  to  indicate  the  texture  of  height  disturbance.  The  texture  of 
bump  represents  a  random  profile  variation  of  height  according  to  a  given  statistical  distribution. 

An  example  of  profile  of  512  x  256  points  is  presented  on  figure  4: 
height  values  could  take  256  values, 
height  RMS  is  38. 

The  ratio  of  15%  (38/256)  will  remain  constant  for  this  texture.  The  peak-to-peak  value  and  the  length  of  the 
motif  (of  512  points)  would  be  the  two  parameters  that  have  to  be  tuned  according  to  a  given  wavelength 
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Figure  5:  Statistica 


profile  parameters 


The  density  of  the  irregularities  of  a  surface  is  described  by  the  autocorrelation  function  that  evaluates  the 
degree  of  covering  between  a  profile  and  its  shifted  image  of  a  distance  R. 

It  is  defined  by  p(R)  =  <7^C( R )  where  C(R)  is  the  function  standardised  autocorrelation  function  or 
coefficient  of  correlation  and  C7  is  the  RMS  height  of  roughness. 

Most  of  profile  modelling  generally  admits  that 

C(R)  =  exp( -( R/ )  for  a  Gaussian  function  and 

C(R)  =  exp(-( R/ L^))  for  an  exponential  function  where  is  the  correlation  length  ,  defined 
as  the  distance  for  which  the  coefficient  of  correlation  is  equal  to  1/e  or  -4.3  dB. 
gives  an  estimate  of  the  size  of  the  irregularities  of  surface.  Correlation  length  of  several  profiles  was 

calculated  and  is  presented  on  previous  figure.  It  is  important  to  remind  that  this  length  of  correlation  does  not 
depend  on  the  direction  chosen  on  the  profile. 

A  standardised  and  centred  histogram  on  256  profiles  of  512  points  extracted  from  this  texture  of  bump  is 
presented  on  the  figure  according  to.  It  shows  the  Gaussian  and  isotropic  behaviour  of  the  distribution  heights 
on  this  texture. 


3.3.  Roughness  effects  on  canonical  objects 

The  roughness  previously  described  is  used  to  measure  its  impact  on  the  RCS  of  simple  canonical  objects 
such  as  plate  planes,  dihedral,  trihedral  and  canonical  scene  (from  SET  069).  The  length  of  the  motif  (of  512 
points)  is  1  m  and  calculations  have  been  carried  out  at  34  GHz.  Typical  results  according  to  the  height  of  the 
texture  of  roughness  are  presented  on  the  figure  6. 
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Figure  6:  Roughness  effect  on  RCS  of  canonical  targets 


On  this  figure,  attenuation  on  the  maximum  RCS  is  observed.  As  an  example,  for  a  peak-to-peak  roughness 
value  of  one  wavelength  X.,  following  values  could  be  noted: 

- 1 5  dB  for  the  plate, 

-8  dB  for  the  dihedral, 

-5  dB  for  the  trihedral. 

Less  attenuation  is  observed  in  the  backscattering  direction  for  the  dihedral  and  the  trihedral  even  though 
multiple  interactions  are  involved.  In  fact,  the  increase  of  the  angle  of  incidence  on  the  constitutive  plates  of 
each  object  decreases  the  effect  of  the  roughness  in  agreement  with  the  modified  Fresnel  coefficient  reminded 
above.  It  is  important  to  note  that  at  the  same  time,  the  sidelobes  level  is  globally  increasing. 


3.4.  Results  on  the  ZSU  23-4  60  kfacets 

On  figure  8,  from  left  to  right  we  respectively  have  a  photo  and  2  pictures  from  the  600  kfacets  and  60  kfacets 
CAD  files.  It  is  important  to  note  that  in  terms  of  EM  interactions  at  Ka  band,  the  2  pictures  do  not  really 
represent  the  same  object,  particularly  concerning  the  caterpillar  tracks. 
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Figure  8:  Representation  of  the  ZSU  23-4 


After  analysis  and  adjustment  of  the  effects  on  the  canonical  objects,  it  was  decided  to  apply  previous  texture 
to  the  60  kfacets  CAD  file  of  the  ZSU  23-4.  First  results  are  presented  below  on  figures  9  and  10. 


Figure  8:  ISAR  images  with  and  without  roughness  texture 


Concerning  the  ISAR  images,  high  returns  from  specular  effects  are  present  on  the  central  image  and  the 
distribution  of  scatterers  is  different.  RCS  investigation  has  been  carried  out  from  these  images. 


Figure  9  presents  the  RCS  behaviour  from  each  target.  Smooth  configurations  of  600  kfacets  and  60  kfacets 
models  are  really  different.  Comparison  between  results  from  smooth  and  rough  60  kfacets  CAD  file  shows 
reduction  of  principal  peaks.  CAD  files.  It  is  more  difficult  to  compare  results  from  rough  60  kfacets  and 
smooth  600  kfacets. 

4.  ONGOING  WORK  -  CONCLUSION 

The  principal  shortcoming  of  a  RCS  calculation  in  MMW  is  the  too  significant  value  of  "specular" 
backscattering  due  to  unsatisfactory  available  CAD  files  compared  to  the  real  targets.  FERMAT  has  the 
potential  to  take  into  account  roughness  texture  thanks  to  the  adaptive  antialiasing.  From  these  first 
experiments  of  roughness  modelling,  it  appeared  the  results  are  satisfying  compared  to  available 
measurements. 

Comparison  on  global  RCS  (frequency  averaged)  and  ISAR  images  on  measurement  references  from  SET  069 
data  pool  and  other  available  measurements  from  DCE/CELAR  will  be  carried  on.  Statistical  analyses  on 
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ISAR  images  evaluated  in  terms  of  peak  values,  mean  values,  quadratic  mean  will  be  assessed  and  could  allow 
an  adjustment,  on  a  given  class  of  vehicle,  of  the  random  part  of  the  scattering  process. 
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fluctuations  through  a  density  of  probability) 

•  Surface  roughness  by  RMS  height  and  eorrelation  length 
characterisation 
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Shooting  and  Bouncing  Rays 


General  Principle 

>  SBR  :  Optimised  technique  to  explore 
geometrical  data  base,  validated  in  optical 
domain 

p-  Set  of  rays  representing  an  incident 
plane  wave 


Shooting 


Shooting  window 


pixel 
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Electromagnetic  Models  :  asymptotic  methods  of  rays 
and  current  on  dielectric  or  metallic  surfaces 


•  Geometrical  Optics  GO:  coupling,  and  forward 
scattering  (multiple  interactions)  ,  very  near  field 

•  Physical  Optics  PO:  scattering  to  observation  (mono 
or  bistatic  configuration) 

•  Uniform  Theory  of  Diffraction  UTD:  (near  field  edge 
diffraction) 

•  Equivalent  Currents  Method  ECM:  (edge  diffraction 
to  observation) 
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Examples  and  Results  on 
Various  Data  Base 
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Radar  Signature  Prediction 


•  Tilted  Trihedral  Comer  (dimensions  -140  cm  x  60  cm  x  20  cm) 

•  Polished  aluminum  full-scale  canonical  target 

•  10  GHz  RCS  data  published 

•  Workshop  on  19  Febmary  2003  in  ONERA  Salon  de  Provence  at 
9,  34,  94  GHz 
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Tilted  trihedron  W-Band  RCS  at  10°  depression; 

Interaetion  decomposition 


WBand-VVRCS 


—  PO 

single 

interaction 

—  GO-PO 

double 

interaction 

—  GO-GO-PO 

triple  interaction 

—  GO-GO-PO 

+  PTD 
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RCS  (dB.m2) 


Tilted  trihedron  W-Band  RCS  at  10°  depression; 

Validation 


WBand-VVRCS 


- FERMAT 

-  XPATCH 
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Radar  Cross  Section  Results  on 
VAB  (Vehicule  de  1’ Avant  Blinde) 
Repair  and  Recovery  Vehicle 
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Radar  cross  section  at  millimetre  wavelength 
Comparison  with  measurements 


CAD  virtual  object 
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Results 


Measurement 

Calculations 
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ISAR  images  at  millimetre  wavelength 


Angular  sector  -0.5/-0.5  deg 

Measurements  Calculations 
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5dB 


Profile  Statistical  Parameters 


Gaussian  distribution  of  the  heights,  characterised  by: 

•  256  levels  (8  bits) 

•  An  extended  area  of  1 

Adjustment  of: 

•  rms  height 

•  correlation  length 

•  antialing  parameters 
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Normalised  height  on  256  levels 


Roughness  motif  on  512  x  256  points 
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Autocorrelation  of  3  bump  profiles  on  512  pts  representative  of  57  wavelengths 


Normalised  distance  relating  to  wavelength 
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3  bump  profiles  of  15%  rms  height,  normalised  by  the  wavelength 
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Normalised  and  centred  histogram  of  heights  on  256  profiles  of  512  points 
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Available  data  and  measurement  on  ground  vehicles  : 

^  Data  pool  group  SET  069 
•ARE 

•  ZSU  CAD  600  kfacets  and  60  kfacets 

•  T72  measurements 

•  ZSU  23-4  FGAN ISAR  measurements 

•  T72  FGAN  measurements 
^  Specific  measurements  on  rough  plates 
^  Others  ... 


First  analysis  of  roughness  effect  on  ground  vehicles  : 

^  Based  on  previous  investigations  on: 

•  a  VAB  (Vehicule  de  I’Avant  Blinde)  Repair  and  Recovery  Vehicle 

•  the  use  of  bump  texture  in  IR  applications 

^  Specific  roughness  analysis  on  canonical  objects  (plates,  dihedral,  trihedral,  test 
scene) 

^  First  applications  on  the  ZSU  CAD  60  kfacets 
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Roughness  effect  on  a  plate  RCS 


0.0  A. 
0.5  A. 
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Roughness  effect  on  a  dihedral  RCS 
peak-to-peak  value  as  parameter 


Angles  in  degrees 
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Roughness  effect  on  a  trihedral  RCS 
peak-to-peak  value  as  parameter 
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Canonical  Target  Imaging  35  GHz  W 


7.6  dBm^ 
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Range  in  m 


3 


Canonical  Target  Imaging  Set  35GHzVV  Lamda/4 
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Range  in  m 
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Canonical  Target  Imaging  Set  35GHzVV  Lamda/4 
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Bump  texture  effect  and  modified  Fresnei 
formulations  (R//^*exp[-2*(  27racos(e)/?i)2] 


Evaluation  tools  (metric  comparison) 


Evaluation  and  comparison  tools  worked  on: 

Statistical  evaluation  of  ISAR  images 

•  peak  values, 

•  mean  value, 

•  quadratic  mean. 

Comparison  on  global  RCS  (frequency  averaged) 
Measurements  references  from  SET  069  data  pool  and 
other  available  measurements  from  DECE/CELAR. 
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ZSU  23-4  600  kfacets 
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* 


ZSU  23-4  600  and  60  kfacets 
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CAOZSU23-4 


CALCULATION 


600  kfacets 


Total  azimuth  angle 


180 


Bandwidth 


MHz 


Step  frequency  number 


256 


Total  angle  number 


12000 


17.25 


Total  caculation  duration 

on  PC 


days 


mn/by  angle 


ISAR  imaging 


Total  ISAR  image  number 


295 
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Meter 


ZSU  23-4  35  GHz  W  dBm^ 
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ZSU  23-4  35  GHz  HV 
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RCS  ZSU  23-4  600k /60k  35  GHz  VV  (frequency  mean) 
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RCSZSU23-4  60k  35  GHz  VV  (frequency  mean) 
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ZSU  23-4  60k  Bump2  35  GHz  VV  12  degres 
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Conclusion  on  the  methodology 


From  the  first  experiment  of  modelling  of  the  VAB  it 
appeared  that  with  the  principal  shortcoming  of  a  RCS 
calculation  in  MMW  is  the  too  significant  value  of  "specular" 
compared  to  a  diffuse  scattering. 

Statistical  analyses  on  ISAR  images  evaluated  in  terms 
of: 


•  peak  values, 

•  mean  value, 

•  quadratic  mean, 

can  allow  an  adjustment,  on  a  given  class  of  vehicle,  of  the 
random  part  of  the  scattering  process. 

Comparison  on  global  RCS  (frequency  averaged) 

Measurements  references  from  SET  069  data  pool  and 
other  available  measurements  from  DCE/CELAR  will  be 
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